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Abstract 

We study the amounts of fine-tuning of the parameters of the MSSM with 
non-universal soft supersymmetry-breaking contributions to the Higgs masses 
(the NUHM) that would be required for the relic neutralino density to lie 
within the range favoured by WMAP and other astrophysical and cosmo- 
logical observations. Such dark matter fine-tuning is analogous to the com- 
monly studied electroweak fine-tuning associated with satisfying the elec- 
troweak symmetry breaking conditions, which we also study for complete- 
ness. We identify several distinct regions of the NUHM parameter space: 
a bulk region, a f — Xi coannihilation region, a pseudoscalar Higgs funnel 
region, a focus-point bino/higgsino region and a P — Xi coannihilation re- 
gion. Within each region, we analyse specific representative points for which 
we provide breakdowns of the contributions to the dark matter fine-tuning 
associated with the different NUHM parameters. In general, the NUHM of- 
fers points with both smaller and larger amounts of dark matter fine-tuning 
than points in the corresponding regions of the CMSSM. Lower amounts 
of dark matter fine-tuning typically arise at points where several different 
(co) annihilation processes contribute, e.g., at junctions between regions with 
different dominant processes. We comment on the prospects for using col- 
lider measurements to estimate the likely dark matter density within the 
NUHM framework. 



1 Introduction 



The primary utilitarian motivation for supersymmetry being accessible to experi- 
ments at the electroweak scale, e.g., at the LHC, depends on its ability to alleviate 
the problem of fine-tuning of electroweak symmetry breaking present in the Stan- 
dard Model [HE]- A supplementary phenomenological motivation for weak-scale 
supersymmetry is its ability to provide the cold dark matter required by astro- 
physics and cosmology P, ]. This is a natural feature of supersymmetric models 
that conserve R parity, with the lightest neutralino x? being particularly well-suited 
to provide the preferred amount of cold dark matter if it is the lightest supersym- 
metric particle (LSP) and weighs less than about 1 TeV [31 [5] . Within the general 
supersymmetric framework, one may find more plausible regions of the supersym- 
metric parameter space that are less fine-tuned, in the sense that the values of the 
model parameters chosen at some high input scale require less delicate adjustment 
in order to obtain the correct value of the electroweak scale [HE], as measured by 
Mz, or the correct value of the cold dark matter density QcDAih"^ [3, El El fTO] . 

It is hard to make this type of plausibility argument at all rigorous: it is notoriously 
difficult to make probabilistic statements about the unique (by definition) Universe 
in which we live, it is largely a matter of personal choice which derived quantity 
one should consider and which input parameters one wishes to avoid fine-tuning, 
it is difficult to argue conclusively for the superiority of one measure of fine-tuning 
over any other, and even less easy to agree on a 'pain threshold' in the amount 
of fine-tuning one is prepared to tolerate [6]. Nevertheless, within a given model 
framework with its specific input parameters, it is legitimate to consider some 
important derived quantity such as QcDMh"^, and compare the amounts of fine- 
tuning required in different regions of its parameter space, which frequently do not 
depend very sensitively on the specific sensitivity measure employed. 

Moreover, even if one does not accept that the less sensitive parameter regions are 
more plausible, measuring the dark matter fine-tuning may have other uses. For 
example, one hopes (expects) some day to discover supersymmetry and start to 
measure the values of its parameters. Unavoidably, these will have non-negligible 
measurement errors, and these uncertainties propagate via the dark matter fine- 
tuning parameters into the calculation, e.g., of QcDMh"^- One of the key features of 
supersymmetry is its ability to provide a calculable amount of cold dark matter, and 
it is interesting to know how accurately which of its parameters must be measured 
in order to calculate flcDnh'^ with an accuracy comparable to that quoted by 
astrophysicists and cosmologists [Z1[II]. An accurate calculation of VtcDuh'^ might 
also reveal some deficiency of the supersymmetric explanation of the cold dark 
matter, and possibly the need for some other new physics in addition. 

It is important to note that the parameters we refer to here are the GUT-scale 
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soft supersymmetry-breaking masses and couplings, whereas experiments would 
measure directly physical masses and mixings at much lower energies. Ideally, one 
would calculate the relic density directly from the low-energy measurements of 
MSSM parameters. However it will be difficult, if not impossible, to pin down all 
the key parameters using LHC data alone, except by making supplementary as- 
sumptions about the pattern of supersymmetry breaking at the GUT scale, as we 
do here. Assuming a structure of GUT-scale unification, one may use experimental 
measurements to constrain these fewer high-energy parameters. The strength of 
the constraints will depend on the magnitudes of these parameters and the ex- 
perimental tools available. Very likely some accelerator beyond the LHC will be 
needed, but we do not yet know what will be available. The fine-tuning measures 
we calculate here show clearly which of the high-energy parameters are most im- 
portant for a precise calculation of the relic density, and hence contribute to the 
'wish list' for such an accelerator. 

For these reasons, we make no further apologies for considering the fine-tuning 
of flcDMh"^ in this paper, which we shall refer to as "dark matter fine-tuning" 
[3, El El [To] to distinguish it from the more commonly studied "electroweak fine- 
tuning" [H [6], which we also consider for completeness. The issue of dark matter 
fine-tuning has been considered previously in the context of several different mod- 
els including the constrained minimal supersymmetric extension of the Standard 
model (CMSSM) [7], in which the soft supersymmetry-breaking scalar masses mo, 
gaugino masses mi/2 and trilinear parameters Aq are each assumed to be uni- 
versal, a more general MSSM with non-universal third family scalars and gaug- 
ino masses [H], a string-inspired non-universal model [H] and SUSY GUTs with 
non-universal gaugino masses [TU]. Here we extend such considerations to models 
with non-universal soft supersymmetry-breaking contributions to the Higgs masses 
(NUHM). Within this NUHM framework, the independent input parameters may 
be taken as [El [131 [H] 

ttNUHM = {mo, niHi, mH2, "^1/2, Aq, tan/3, sign(/i)} , (1) 

and we take as our measure of dark matter fine-tuning the quantity 



An = Maxi 



fi^ dai 



(2) 



Our objective will be three-fold: to compare the amount of dark matter fine-tuning 
required within the NUHM to that required within the CMSSM, to identify the 
regions of the NUHM parameter space that require relatively less (or more) dark 
matter fine-tuning, and thereby to quantify the accuracy in the determination of 
the GUT-scale NUHM parameters that would be needed in order to calculate VL^h? 
with any desired accuracy. 



2 



The regions of the NUHM parameter space where Q^h"^ falls within the range 
favoured by WMAP and other experiments has been studied quite extensively, for 
example in [H]. It shares several features in common with the more restrictive 
CMSSM framework proposed in [15] and extensively studied in [16]. For example, 
there are regions where Xi - stau coannihilation is important, and others where Xi 
pairs annihilate rapidly via direct-channel H, A poles. However, other possibilities 
also occur. For example, there are regions where Xi ' sneutrino coannihilation is 
dominant. Also there are regions where rapid-annihilation and bulk regions, which 
are normally separated by a coannihilation strip, approach each other and may even 
merge. As we discuss below in more detail, the sneutrino coannihilation regions 
exhibit relatively high dark matter fine-tuning, whereas the 'merger' regions may 
require significantly less dark matter fine-tuning. 

In this work we provide a first calculation of the dark matter fine-tuning for the 
regions of the NUHM that are favoured by dark matter measurements. In addition, 
we present a first calculation of the electroweak fine-tuning within this model and 
update the parameter scans for the current measurement of the top mass. 

The rest of the paper is laid out as follows. In Section [2] we summarise the methods 
used in our numerical studies. Next, in Section [3] we review the familiar case of 
the CMSSM, which serves as a baseline for later comparison. Then, in Section H] 
we study dark matter within the NUHM model in which universality between the 
soft supersymmetry-breaking masses of the sfermions (squarks and sleptons) and 
Higgs multiplets is broken. Finally, in Section [3] we present our conclusions. 

2 Methodology 

2.1 Codes 

In order to study the low-energy phenomenology of the NUHM, we need a tool 
to run the mass spectrum from the GUT scale down to the electroweak scale 
using the renormalisation group equations (RGEs)[T7]. For this purpose we use 
the ROE code SoftSusy [18]. This interfaces with the MSSM package within 
micrOMEGAs [T9], which we use to calculate the dark matter relic density QcDMh"^, 
BR{h — > 57) and 5a^. We take rrit = 170.9 GeV throughout. 

2.2 Theoretical, Experimental and Cosmological Bounds 

After running the mass spectrum of any chosen model parameter set from the 
GUT scale down to the electroweak scale, we perform a number of checks on the 
phenomenological acceptability of the point chosen. A point is ruled out if: 
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1. It does not provide radiative electroweak symmetry breaking (REWSB). Such 
regions are displayed in light red in the subsequent figures. 

2. It violates the bounds on particle masses provided by the Tevatron and 
LEP 2. Such regions are displayed in light blueQ. 

3. It results in a lightest supersymmetric particle (LSP) that is not the lightest 
neutralino. We colour these regions light green. 

In the remaining parameter space we display the 1- and 2-0" regions for 5a^ and 
BR{h — > 57), as well as plotting the 2-o" region for the relic density allowed by 
WMAP and other observations. 

2.2.1 5a^ 

Present measurements of the anomalous magnetic moment of the muon deviate 
from theoretical calculations of the SM contribution based on low-energy e^e~ 
data0. Taking the current experimental world average and the state-of-the-art SM 
value from [21], there is a discrepancy: 

{a^)e.p - {a^)sM = Sa^ = (2.95 ± 0.88) x 10~', (3) 

which amounts to a 3.4-0" deviation from the SM value. As already mentioned, we 
use micrOMEGAs to calculate the SUSY contribution to {g — 2)^. The dominant 
theoretical errors in this calculation are in the SM contribution, so in our analysis 
we neglect the theoretical error in the calculation of the SUSY contribution. 

2.2.2 BR{b 57) 

The variation of BR{b — >■ 57) from the value predicted by the Standard Model 
is highly sensitive to SUSY contributions arising from charged Higgs-top loops 
and chargino-stop loops. To date no deviation from the Standard Model has been 
detected. We take the current world average from [22], based on the BELLE [23] . 
CLEO [21] and BaBar [25] measurements: 

57) = (3.55 ± 0.26) X 10"^ (4) 

^Thc current LEP 2 bound on the lightest MSSM Higgs stands at 114.4 GeV. However, there 
is a theoretical uncertainty of some 3 GeV in the determination of the mass of the light Higgs [50] . 
Rather than placing a hard cut on the parameter space for the Higgs mass, instead we plot a line 
at rrih = 111 GeV and colour the region in which ruh < 111 GeV in very light grey-blue. 

^ There is a long-running debate whether the calculation of the hadronic vacuum polarisation 
in the Standard Model should be done with e~^e~ data, or with r decay data. The weight of 
evidence indicates the e~^e~ estimate is more reliable so we use the e^e~ value in our work. 
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Again, we use micrOMEGAs to calculate both the SM value of BR{b —>■ 57) and 
the SUSY contributions. It is hard to estimate the theoretical uncertainty in the 
calculation of the SUSY contributions, but note that there is an uncertainty of 10% 
in the NLO SM prediction of BR{h — ^ 57) j26] B As with 5a^, we plot the l-a and 
2-cr experimental ranges, and do not include a theoretical error in the calculation. 

2.2.3 VtcDhih'^ 

Evidence from the cosmic microwave background, the rotation curves of galaxies 
and other astrophysical data point to a large amount of cold non-baryonic dark 
matter in the universe. The present measurements [28] indicate the following value 
for the current cold dark matter density: 

f]^^^^/i2 = 0.106 ±0.008. (5) 

We calculate the relic dark matter density with micrOMEGAs using the fast ap- 
proximation. Given a low-energy mass spectrum, micrOMEGAs gives an estimated 
precision of 1% in the theoretical prediction of the relic density. This is negligible 
compared to the present observational error, so the 2-a band plotted takes into 
account only the experimental error B 

In the following Sections, we calculate the dark-matter fine-tuning for any point 
that lies within the 2a allowed region, and indicate the amount using colour coding. 
We also display electroweak fine-tuning contours over the different regions. 

^We recall that micrOMEGAs calculates the SM contribution to BR[b sj) to NLO. A first 
estimate of the SM prediction of BR{b — *■ 57) to NNLO was presented in [57] ■ This showed a 
decrease of around 0.4 x 10^^ in the central value of the SM prediction. The implementation of 
the NNLO contributions in the calculation is non-trivial and its implementation in micrOMEGAs 
is currently underway. As a result we do not include this decrease in the results we present, 
but instead note that positive SUSY contributions to BR{b 57) look likely to be favoured in 
future. This would favour a negative sign of fi and thus cause tension with (g — 2)^. 

"^Wc emphasize that the quoted 1% accuracy is for a given low-energy spectrum, which is 
obtained using sof tsusy. However, there are differences in the details of the mass spectrum 
between codes |29| . for given high-energy inputs, and different dark matter regions have different 
levels of sensitivity to these variations: see [3D] for a detailed study. The result of the discrepancies 
between codes is to move the dark matter regions slightly in the GUT scale parameter space. 
As we are interested in broad features of these regions, rather than their precise locations, these 
uncertainties are not important for our purposes. 
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3 The Constrained Minimal Supersymmetric Stan- 
dard Model 



We first review the familiar Constrained Minimal Supersymmetric Standard Model 
(CMSSM) [15], [16]. wfiicfi serves as a standard to wliicli we compare tlie parameter 
space of the NUHM. 

The CMSSM has a much simpler spectrum of soft masses than the full MSSM. 
First, all of the soft squark and slepton (mass)^ matrices are chosen to be diagonal 
and universal at the GUT scale with the diagonal entries equal to mg. Secondly 
both the soft Higgs (mass)^ are also set equal to ttIq. Additionally, all the gaugino 
masses are assumed to be unified with a value mi/2 at the GUT scale. Finally, 
we take the trilinear coupling matrices to have only one non-zero entry (the third- 
family dominance approximation) and assume that all these entries are equal to a 
common value Aq. Requiring that electroweak symmetry be broken radiatively to 
give the observed electroweak boson masses, we trade the soft parameters fi and 
B for tan f3, the ratio of the Higgs vevs, and the sign of /x. This results in a model 
with four free parameters and a sign: 

acMSSM & {mo, mi/2, Aq, tan/?, sign(/i)}. (6) 

Although our main focus is the dark-matter fine-tuning, we also report the required 
amounts of electroweak fine-tuning for specific cases of interest. 

In Fig. [1] we show the {1^0,7711/2) plane of the CMSSM for Ao = 0, tan/? = 10 and 
sign(/x) positive. At low mo the parameter space is ruled out because < m^o, 
resulting in a stau LSP (light green). Regions at low mi/2 are ruled out by LEP 2 
bounds on the masses of the charginos and sleptons (light blue). A larger range 
of mo is ruled out by the absence of a light Higgs boson, and we shade the region 
with mh. < 111 GeV (light grey with a black boundary). Finally, the model fits 
the current measurement for low mo and mi/2- The la and 2a bounds are 
shown as a short and long-dashed green lines respectively. At larger mo and mi/2, 
the SUSY contribution becomes small as the sparticles that contribute in loops 
become heavy. Therefore 6a^ ^ favours relatively light soft masses. 

In the remainder of the parameter space we find the regions that fit the WMAP 
strip at 2a, and for each such point we calculate the fine-tuning of the dark matter 
density. Each point is then plotted with a colour that corresponds to the value of 
the tuning via the log-scale on the right hand side. 

There are two distinct regions that fit the WMAP measurement of the relic density. 
The first region lies at mi/2 ~ 200 GeV, mo ~ 70 GeV and contains the point Al. 
This is the bulk region, which is adjacent to the light blue region that is ruled 
out by the LEP 2 slepton mass constraint. The lightness of the sleptons enhances 
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m\/2{Ge\'} 

Figure 1: The (mi/2;"io) plane of the CMSSM with Aq = 0, tan/? = 10 and sign(/i) 
+ve. 

neutralino annihilation via t-channel slepton exchange to the extent that it allows 
bino dark matter to fit the WMAP relic density measurement. This process is 
relatively insensitive to the precise masses of the neutralino and the sleptons. This 
is reflected in the fact that much of the bulk region is plotted in yellow, signifying 
< 1. From the breakdown of the tunings of point Al in Table [T] it is clear 
that the tuning is mainly in the parameters mi/2 and rriQ. The sensitivity to mo is 
to be expected, as t-channel slepton exchange is clearly dependent on the mass of 
the exchanged particle. This also explains the sensitivity to mi/2 since, although 
the masses of the sleptons are determined by mo at the GUT scale, the running 
to low energies is dominated by mi/2- Thus, in this region the masses of the light 
sleptons are sensitive to both mg and mi/2- The relic density is also dependent 
upon the mass of the LSP, in this case (mainly) a bino with a mass determined 
primarily by mi/2- It is therefore apparent why the sensitivity of the bulk regions 
lies primarily with mo and mi/2. The fact that no precise cancellations or balancing 



7 



Parameter 


Al 


A2 


value 




value 




mo 


60 


0.62 


100 


5.7 


mi/2 


200 


0.99 


500 


5.8 


tan (3 


10 


0.13 


10 


1.5 


An 




0.99 




5.8 


Aew 




37 




190 



Table 1: A summary of the properties of points Al and A2, shown in Fig. [H chosen 
as representatives of the bulk region (Al) and stau-coannihilation region (A2). We also 
present a breakdown of the dark-matter and electroweak fine-tunings with respect to 
each parameter of the CMSSM, except for Aq, which we fix: Aq = here and elsewhere 
in this paper. 

of parameters is required explains why the tuning is low. However, the bulk region 
lies inside a light grey region signifying a Higgs with a mass less than 111 GeV. 
Therefore we do not consider this region further here. 

The second region that fits the observed relic density is the band that contains 
point A2. This band lies alongside the area that is ruled out by a stau LSP (light 
green). Along the edge of the light green region, the mass of the stau is close to 
that of the lightest neutralino, resulting in comparable number densities of both 
particles around the time of freeze-out. Thus many more annihilation channels 
must be considered, such as the annihilations of f — f and f — Xi in addition to 
Xi~Xi- This suppresses the number density of SUSY particles and thus reduces the 
resultant dark matter relic density. The effect of coannihilation depends strongly 
upon the number density of NLSPs at freeze-out. This is in turn very sensitive 
to the mass difference between the NLSP and the LSP. Thus, it may at first sight 
seem surprising that the f coannihilation band is plotted here in green and red, 
corresponding to relatively low dark matter fine-tuning = 3 — 11. The reason, 
as discussed in [H], is that in this region the RGE running results in the mass of 
the right handed f and the mass of the lightest neutralino both being dominantly 
dependent on mi/2, with the stau having a smaller secondary dependence on ttiq. 
As a result, their masses vary together as the soft masses are varied, and Am 
is remarkably insensitive to mo or mi/2. This mitigates the normal sensitivity of 
the coannihilation region. We take the point A2 as a representative point in this 
band and provide a breakdown of the individual tunings in Table [TJ The tuning is 
equally dependent upon mo and mi/2, as expected. 

As well as the dark-matter fine-tuning, for reference we also calculate the fine- 
tuning required to fit the electroweak boson masses. This calculation is performed 
across the parameter space using the same measure as we use for the dark matter 
fine-tuning. We plot the tuning using blue dot-dashed contours in Fig. [1], and label 
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each contour with their respective values of ISP"^ = max(Af^). We also list the 
value of in the last row of Table [T] for both points we consider. The W boson 
mass is in general the result of a careful balancing act between the soft masses: 
the larger the soft masses, the more precisely they must cancel to give the W 
boson mass. Therefore it is unsurprising that the amount of fine-tuning required 
to obtain the correct electroweak symmetry breaking increases smoothly with the 
increasing soft masses. The electroweak and dark matter fine-tunings are largely 
independent. 




miy2(GeV) 

Figure 2: The (mo, rni/2) plane of the CMSSM parameter space for ^0 = 0, tan/3 = 50 
and sign(/i) +ve. 

Other regions in which the CMSSM fits the observed relic density are seen when 
we consider larger values of tan /3. In Fig. [2] we consider a wider range in the 
(mo,mi/2) plane and take tan/3 = 50. 

Low mo is once again ruled out because the stau is the LSP (light green) . Equally 
low mi/2 and mo are ruled out by searches for the lightest Higgs (light grey with 
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Parameter 


A3 


A4 


A5 


A6 


value 




value 




value 




value 




mo 


2030 


27 


540 


5.0 


277 


19 


1400 


13 


mi/2 


500 


18 


600 


8.0 


350 


12 


250 


22 


tan f3 


50 


2.0 


50 


76 


50 


48 


50 


8.2 






27 




76 




48 




22 






150 




230 




92 




48 



Table 2: Here we summarize properties of the points A3-6, shown in Fig. [21 which 
represent the higgsino/bino region (A3, 6), the pseudoscalar Higgs funnel (A4) and the 
stau-coannihilation region at large tan/3 (A5). We also present a breakdown of the 
dark matter fine-tuning with respect to each parameter, and also list the electroweak 
fine-tuning. 

a black boundary) or LEP 2 bounds on the chargino and light sleptons (solid light 
blue). Finally large mo and low mi/2 are ruled out as they do not give REWSB 
(light red). We plot the different regions that fit the relic density using the same 
colour coding for dark-matter fine-tuning as was used in Fig. [T] 

We now distinguish four regions in which the CMSSM can account for the observed 
relic density of dark matter. The least interesting of these regions is the bulk region, 
which once again appears at low mo and low mi/2, but requires a Higgs lighter than 
current search limits allow. Moving on from the inaccessible bulk region, we again 
find the coannihilation strip along the side of the region ruled out because the 
stau is the LSP. In comparison to FigH], it is thin and broken. This is a sampling 
artefact, arising because we have extended the range of the mo scan, and the 
resolution is not fine enough to resolve the full strip. In contrast to Fig{Tl here 
the strip is plotted in purple, designating a dark matter fine-tuning in the range 
= 30 - 100. 

Consider now the representative point A5 and the corresponding fine-tuning. At 
large tan (3, the mass of the f^j runs down more than for low tan f], and therefore a 
larger soft mass is necessary to avoid a f LSP. For low tan (] the f coannihilation 
region occurred when mo -C mi/2, whereas here mo and mi/2 are much closer in 
magnitude. This reduces the dominance of mi/2 in the running and restores the 
need for a precise balance of mo and mi/2 to keep the r — x? mass difference small. 
However, it is the influence of tan (3 over the running that dominates the sensitivity. 

The second region of interest shows up as two diagonal bands starting beside the 
light green region ruled out because the stau is the LSP. These lines lie to either 
side of the line along which 2m^o = m^, and neutralino annihilation proceeds via 
the resonant production of a pseudoscalar Higgs boson. This region is known as 
the pseudoscalar Higgs funnel region. The WMAP lines lie along either side of this 
resonance, where there is just enough of an enhancement in the annihilation cross 
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section to allow a bino LSP to account for the observed relic density. As a result, 
the dark matter fine-tuning price of such a region is large. This is refiected in 
the purple and blue shading of the funnel region, corresponding to a dark matter 
fine-tuning in the range = 30 — 300. 

We consider the representative point A4 and break down the fine-tuning. The 
mass of the pseudoscalar Higgs is dependent upon the details of the running of 
the soft Higgs mass-squared terms, which also determine the Higgs vevs and thus 
tan/3. As we require REWSB and set tan/3 at the start, we can run this chain 
of logic back the other way. A fixed value of tan j3 requires a specific value of the 
soft Higgs mass-squared terms at the low-energy scale. Thus the value of tan/3 
has a significant impact on the pseudoscalar Higgs boson mass. This explains the 
dominant sensitivity of the pseudoscalar Higgs funnel region to tan/3. This is to 
be contrasted with the NUHM, which will allow access to the pseudoscalar Higgs 
funnel region without large tan /3, so that naively one might expect a more natural 
funnel region. 

The final WMAP strip lies at large mg in all panels. The corner of the parameter 
space at large mo and low mi/2 (shown in light red) is ruled out as /i^ < 0, signifying 
a failure of REWSB. Along the perimeter of this region ~ 0, and decreases 
steadily as one approaches the boundary. As /i decreases, the higgsino component 
of the LSP increases. Higgsino dark matter annihilates very efficiently resulting 
in VtcDKih^ ^ ^cdm'^^^' whereas bino dark matter generally gives flcoMh^ ^ 
^ODM^^^- Along the WMAP strip, the higgsino component of the Xi is large 
enough that the annihilation cross section is enhanced sufficiently for the relic 
density to fit the WMAP data. This region is known as the focus-point region. It 
is sensitive to the composition of the lightest neutralino, and so depends upon the 
difference between /i and Mi at the electroweak scale. The region is plotted in red 
and purple, with a couple of green regions at low mi/2- This corresponds to a dark 
matter fine-tuning of «i 10 at low mi/2, rising to a dark matter fine-tuning of 
A*^ 60 at large mo and mi/2- The kink at mi/2 ~ 450 GeV, where the tuning 
drops (signified by a change from purple to red) corresponds to the top quark mass 
threshold where processes of the form XiXi ~^ tt become kinematically allowed. 

We take points A3 and A6 as representative points in this region and break the 
dark matter fine-tuning down into its individual elements. At point A3 the LSP 
is predominantly bino with a small, but significant, higgsino component. The hig- 
gsino component is determined by the relative size of fi{EW) to Mi{EW). As we 
set fi{EW) by requiring REWSB, fi{EW) is determined by the running of the soft 
Higgs masses. Thus fi{EW) is sensitive to the soft Higgs mass-squared terms at the 
GUT scale (set to mg in the CMSSM). It is also sensitive to strongly-interacting 
sparticle masses through the RGEs. This once again brings in a sensitivity to 
mo, but also to M3 through its strong infiuence on the stop quark mass. There- 
fore fi{EW) is sensitive to both mo and mi/2- To achieve the correct balance of 
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bino and higgsino components one needs to balance fi{EW) and Mi{EW). These 
both depend strongly upon mi/2, with /i also dependent upon rriQ. This common 
dependence on mi/2 reduces the dark matter fine-tuning below what one would 
expect. 

At point A6 we consider the kink at the bottom of the higgsino/bino line, where 
the higgsino and bino components are almost equal. This would normally re- 
sult in extremely efficient annihilation of neutralinos in the early universe and 
give VLcDuh^ VL'^^^j^^h? . However at this point the neutralino mass is m^o = 
79.6 GeV, disallowing the annihilation of neutralinos to Higgs, Z bosons and ti. 
It is also on the W boson threshold, suppressing annihilation to W^. With these 
annihilation channels ruled out, the normally efficient annihilation of an LSP with 
a substantial higgsino fraction is suppressed sufficiently to fit the WMAP relic den- 
sity. As with point A3, the determining factors for the annihilation cross section 
are the relative sizes of ^{EW) and Mi{EW), resulting in a similar pattern and 
magnitude of dark matter fine-tuning measures. 

Once again we also calculate the electroweak fine-tuning across the plane and 
plot the corresponding contours of A^^. As in the case of Fig. [H the degree of 
electroweak tuning increases steadily with increasing mi/2- This is because the 
dominant soft term is M3, the gluino mass. As mi/2 increases, the gluino mass and 
squark masses increase, requiring more precise cancellations to reproduce the W 
mass, and thus greater electroweak fine-tuning. 



Region 


Tuning Range 


Pseudoscalar Higgs Funnel 
f — x\ coannihilation (large mg, ?77,i/2 or tan (3) 
f — x\ coannihilation (low mo, mi/2 and tan/5) 
bino/higgsino region 


60-1200+ 
30-60 
3-10 
10-60 



Table 3: Here we summarize the DM annihilation channels present within the CMSSM 
and their associated tunings. This provides a reference to which we will compare the 
regions accessible within the NUHM. 

Now that we have considered the different dark matter regions present within 
the CMSSM, we summarize the typical tunings in each case. We distinguish in 
Table E] the principal dark matter regions present within the CMSSM and the 
corresponding amounts of dark matter ffne-tuning. With this as our starting point, 
we now consider the case of the NUHM model in which the universality between 
the soft Higgs masses and the soft sfermion masses is broken. 
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4 The NUHM 



We now consider the MSSM with non-universal Higgs soft masses (NUHM) [121 [131 
m]. After breaking the universahty between the soft Higgs and sfermion masses of 
the CMSSM we have the following independent inputs: 

aNUHM = {mo, rriH^, tuh^, mi/2, Aq, tan/5, sign(/i)}, (7) 

where ttlhi and characterize the independent soft Higgs masses. These are 
subject to constraints arising from vacuum stability and cosmological considera- 
tions, and may be negative. As long as m|^^ ^ + fi^ > at the GUT scale, there 
is no dangerous high-scale vacuum state, but specifying the precise boundaries of 
the NUHM parameter space lies beyond the scope of this work. 

As with the CMSSM, the NUHM contains a finite number of distinct regions in 
which it can provide the observed dark matter relic density, which were catalogued 
in [ll]. Here we follow the approach of this previous work and reproduce the 
same regions of the parameter space. The plots we present here show the updated 
parameter space for the current world average for the top mass, rrit = 170.9 GeV, 
and include the current dark matter and 6a^ constraints. However, the primary 
goal of this work is rather to analyse the fine-tuning of the dark-matter regions 
of the NUHM. To this end we calculate and plot the dark-matter fine-tuning in 
the allowed parameter space, and also make some comments on the amount of 
electroweak fine-tuning. 

As the NUHM contains the CMSSM as a limiting case, all the dark-matter regions 
present in the CMSSM are present in the NUHM. In addition, there are four new 
regions that are not present in the CMSSM: 

• A pseudoscalar Higgs funnel at low tan/?. 

• A bulk region where Xi annihilation is dominantly mediated via t-channel f 
exchange which does not violate Higgs mass bounds. 

• A i) — Xi coannihilation region. 

• A mixed bino/higgsino region at low mo. 

We shall be particularly interested in understanding how finely tuned the NUHM 
parameters must be in each of these new regions. 
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4.1 Comparison with the CMSSM 



The NUHM contains all the points in the CMSSM parameter space. Therefore, 
we start by studying the tuning of the dark matter points Al-6, presented in 
Tables [Tp| with respect to the parameters of the NUHM. 



Parameter 


Al 


A2 


A3 


A4 




value 




value 




value 




value 




mo 


60 


0.62 


100 


5.7 


2030 


200 




540 


8.1 




60 


0.017 


100 


0.26 


2030 


14 




540 


28 




60 


0.014 


100 


0.26 


2030 


230 




540 


30 


mi/2 


200 


0.99 


500 


5.8 


500 


18 




600 


8.0 


tan P 


10 


0.13 


10 


1.5 


50 


2.0 




50 


76 


An 




0.99 




5.8 




230 




76 






37 




190 




1300 




230 






Parameter 


A5 


A6 
















value 




value 














mo 




277 


23 


1400 


230 












m^i 




277 


1.5 


1400 


10 
















277 


2.5 


1400 


73 












mi/2 


350 


12 


250 


22 












tan P 


50 


48 


50 


8.2 












An 




48 




230 












Aew 




92 




600 









Table 4: A re-analysis of the representative points Al-6 from Figs. \1\2\ calculating their 
tunings with respect the NUHM rather than the CMSSM. 

We show the dark matter fine-tuning of these points with respect to the parameters 
o-NUHM in Table m Point Al represents the bulk region of the CMSSM, which is 
inaccessible because the Higgs is light. The primary annihilation channel is t- 
channel slepton exchange, and the sensitivity in the CMSSM is primarily due to 
mo and mi/2 as they determine the neutralino and slepton masses. This is also 
true in the NUHM, and the sensitivities to the Higgs soft masses are negligible. 

Point A2 represents the low-tan /3 coannihilation region of the CMSSM, in which 
the primary sensitivities were to mi/2 and mo, as these determine the stau mass 
and the neutralino mass. Once again, this picture changes very little in the NUHM, 
with the sensitivity to the soft Higgs masses being negligible. 

Points A3-6 have large tan (3. We recall that A3 and A6 lie in the higgsino-bino 
focus-point region. In the CMSSM the primary sensitivities were to mg and mi/2, 
as mo (and to a lesser extent M3) determine the size of /i, and mi/2 determines 
Mi{EW). Therefore these two parameters determine the mass and composition 
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of the lightest neutrahno, and the total CMSSM dark matter fine-tuning of the 
point in the CMSSM was = 27. In the NUHM we have a very different 
picture. Here the total dark matter fine-tuning is = 230, and the primary 
sensitivities are to mo and ■ This can be explained by the process of radiative 
electroweak symmetry breaking. For electroweak symmetry breaking to occur, 
the Higgs (mass)^ must become negative. By requiring this process to give the 
correct electroweak boson masses we set the size of fi, and thus the magnitude 
of the higgsino component of the lightest neutralino. Therefore to understand 
the sensitivity of a higgsino-bino dark matter region, we must look for the terms 
that contribute to the running of the Higgs mass-squared. First there is the soft 
Higgs mass at the GUT scale, and then there are the running effects, primarily 
the contribution from the stop mass. In the CMSSM, these two terms are coupled, 
reducing the dependence on either one individually. Therefore even though the 
scalar masses are large, the sensitivity of fi to mo remains small. In the NUHM 
there is no connection between the soft sfermion masses and the soft Higgs masses, 
therefore the sensitivity returns. Therefore one should not expect natural bino- 
higgsino dark matter at large mg in the NUHM. The significant increase in the 
electroweak fine-tuning for these points is due to exactly the same physics. 

Points A4 and A5 represent the pseudoscalar Higgs funnel and the stau-coannihilation 
band. At this value of tan /5, the primary sensitivity is to tan (3, a feature not al- 
tered by breaking the universality amongst the scalars. 



4.2 Detour: RGE behaviour with negative masses-squared 

To understand the dependence of the dark matter phenomenology on the NUHM 
GUT scale parameters we need to understand how the soft Higgs masses affect the 
RGEs, and through them the low-energy parameters. Four low-energy parameters 
in particular are useful to consider when we talk about dark matter: fi, m^, and 
tl^r. The higgsino component of the LSP is determined by /i, determines 
the position of the pseudoscalar Higgs funnel, and the lightest stau (a mixture of 
tl^r) mediates the prevalent t-channel slepton exchange annihilation diagrams and 
determines the efficiency of f coannihilation channels. 

After EW symmetry breaking we can write /i as: 

mjr — mjr tan^ (3 1 



^ tan2 p-1 2 ^ ^ ' 

Clearly fi depends on the soft Higgs mass-squared terms and tan P, as well as other 
soft parameters through the RGEs. It is also useful to consider the limit of large 
tan f] where we can approximate ([H]) as: 
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assuming I^tij^jjI ^ ^'z- Therefore for large tan/3, to achieve REWSB and have 
/i^ > we require either negative m^^^ , or very large positive m|^^ . 

The pseudoscalar Higgs mass is determined after EWSB by the relation: 

m\=m^H^+m]j^+2^?. (10) 

Clearly m\ is strongly dependent upon the soft Higgs mass-squared terms, tan/3 
through its effect on /i, and other soft terms through their influence on the Higgs 
RGEs. 

We now consider the explicit form of the soft Higgs mass-squared RGEs: 

= ^A-^9lMl-glMl + hl{ml^+ml^ + ml^+Al) 



+2>hi{ml + + m^^ + A^) - 25j , (11) 
= ^ [-^Ml - glMl + ?,hl{ml + m?^ + m% + A^) + 2^) (12) 



H2) 



where S is definedly: 



2 



fl'l/2 2 , n I 2 2 

b = — [TnTT^—mrT + 2 — m 



2^L + "^Jh + ^e-«) 



-I- (m^ -m| -2m? + m? + m? I). (13) 

The only parameters in these RGEs that we are not free to set at the GUT scale 
are the Yukawa couplings hi. These are set by the requirement that the Higgs vevs 
should give the correct SM particle masses: 

mr,b = -^hr,bVi, mt = "^^t^2- (14) 

Therefore tan/3 influences the RGEs indirectly through its determination of the 
size of the Yukawa couplings. The Yukawa couplings multiply the contribution to 
the RGEs from the soft squark and slepton mass-squared terms and the soft Higgs 
terms. Therefore varying the Yukawa couplings has a large impact on the running. 
As we increase tan (3, we increase V2 with respect to vi, and so we must decrease 
ht and increase hr^b- Therefore we reduce the Yukawa contribution to the running 
of fnjj^, while increasing the contribution to the running of m^^. 

Now consider the RGEs for the right and left handed stau masses: 



d{mlj _ 1 



dt 87r2 



3g^,M^ - glMl + hi (m|^^ + m^^ + m^^ + A^) - 25^15) 
-AglMl + 2hl (m|^^ + m^^ + m^^ + A^) + 4^) (16) 



dt 87r2 

^t'r) — ^ I A^^n/r"^ I o;,2 1^2 , ^2 , ^2 , a2 
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In both cases m'j^^ provides a substantial contribution to the running, with a co- 
efficient of hr- As we have seen, increasing tan/5 increases hr and thus increases 
the impact of the Higgs masses on the running of the staus. Therefore we expect 
any effects of non-universal soft Higgs masses on the stau running to be amplified 
for large tan/3. In the CMSSM, m|^^ will remain positive from the GUT scale 
to the EW scale. Indeed, it is harder to achieve REWSB if m|^^ runs negative. 
Therefore generally this term provides a positive contribution to both the left and 
right handed stau RGE and acts to suppress the stau masses. 

In the CMSSM this poses a problem. As we increase tan/5 we must increase the 
soft stau mass to avoid it becoming the LSP. However as we increase mo we are also 
increasing m'jj_^, and thus increase the effect on the running. This can be avoided 
in the NUHM. We can set small and so avoid a very light f. 

However, there is another more subtle effect. The interaction of the neutralinos 
with the stau also depends upon the composition of the lightest stau which is 
determined by the mixing between tl^r- This mixing is increased if the two states 
are close in mass. In the CMSSM 5* is negligible and so d{mf )/dt^ d{m'j )/dt, 
resulting in the right handed stau always being considerably lighter than the left- 
handed stau. In the NUHM we can avoid this by having a large negative S. This 
acts to suppress the left handed stau mass while increasing the right handed stau 
mass. As we increase the component of the left-handed stau, we increase the 
annihilation rate of neutralinos via t-channel stau exchange. 

4.3 Sample (mo,mi/2) planes in the NUHM 

Having analysed the CMSSM points from the perspective of the NUHM, we now 
turn to a sampling of the full NUHM parameter space. In Fig. [31 we show {mi/2, ttio) 
planes for tan/5 = 10, Aq = and sign(/i) positive. We set the electroweak scale 
parameters /i and to different discrete values in each panel as explained in the 
figure caption. 

As we saw in the previous section, fi and niA are not high-scale inputs into the 
theory, rather they are the low-energy numbers derived from a given set of the true 
input parameters. However, displaying results as functions of these parameters 
can be more informative. As both have a strong dependence on m|^^ ^, we can fit 
a particular value of fi, ijia with the correct choice of mf^^ ^ at the GUT scale. 
Therefore we use a code that varies m|^^ ^ across the parameter space to fit the 
designated low-energy values of fi and rriA- All fine-tunings are calculated in terms 
of the inputs of the NUHM as listed in ([7]). 

By starting with (mo,mi/2) planes, we make contact with the parameter space of 
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300 wn 9at 12m 1500 300 goo 900 120a ism 

mi,j(GeV) i7T,„(GeV) 



Figure 3: The (mo, '^■1/2) plane of the NUHM parameter space with Aq = 0, tan/? = 10 
and sign(^) +ve. The values of fi and vary between the panels: (a) /i = 400 GeV, 
ruA = 400 GeV, (b) = 400 GeV, tua = 700 GeV, (c) /i = 700 GeV, tua = 400 GeV, (d) 
/i = 700 GeV, niA = 700 GeV. This figure can be compared directly to Fig. 2 in [H]. The 
Roman cross in panel (b) indicates the single point where the parameter space makes 
contact with the CMSSM. 

the CMSSM as displayed in Figs. [H [2]0. As before, low mo is ruled out by a f LSP 
(light green), and low mi/2 results in a Higgs with rrih < 111 GeV (light grey with 
black boundary). As before, (5a^ favours light sleptons, and thus low mo and mi/2- 

The dark matter phenomenology shows some similarities to and some marked dif- 
ferences from the CMSSM. First, we see a familiar f coannihilation region alongside 
the region with a f LSP. As in the CMSSM, this region is plotted in red and green, 
designating a tuning of = 3 — 30. The only new feature of the coannihilation 
region here is that effects of the non-universal Higgs soft masses alter the running 
of the stau mass, which allows access to regions with mo = 0. We can access 
mo = with small negative m|^^ and larger negative m'jj^. The combination of a 
small Yukawa contribution (due to low tan/3 along with small |m^J) along with 

^ We note in panel (b) of Fig. [3] the appearance of a CMSSM point, the only point in any of 
these planes where full GUT-scale universality is recovered. 
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J. CLl djlilC u C^l 


Bl 


B2 


B3 


B4 




valup 




valiip 




valup 




valiip 




Tn,r\ 


500 


32 


500 


8.6 


500 


4.6 


500 


12 


2 


-80249 


16 


-126930 


12 


-248480 


0.61 


90625 


2.3 


m?Tr 


461380 


62 


675760 


25 


1202900 


24 


1194100 


60 


mi/2 


435 


39 


540 


19 


750 


18 


750 


38 


tan/3 


10 


5 


10 


3.2 


10 


1.1 


10 


2.9 


An 




62 




25 




24 




60 






150 




220 




390 




390 


/i 


400 




400 




400 




400 






400 




400 




400 




700 





Paramptpr 


B5 


B6 


B7 




value 




value 




value 




mo 


500 


40 


100 


6.3 










-416350 


110 


-400510 


12 


-79656 


2.2 




-24320 


4.1 


-332200 


10 


-266010 


7.4 


mi/2 


442 


52 


400 


3.5 


445 


4.3 


tan j3 


10 


5.8 


10 


0.55 


10 


1.9 


An 




110 




12 




7.4 


A-EW 




250 




250 




250 




700 




700 




700 




rriA 


400 




400 




700 





Table 5: Analysis of the points Bl-7, shown in Fig. O which are representative of 
the pseudoscalar Higgs funnel (Bl,2,4,5), mixed bino-higgsino dark matter (B3) and f 
coannihilation regions (B6,7). We present a breakdown of the dark matter fine-tuning 
with respect to each parameter of the NUHM. We give the value of w?^^ ^ , but the tunings 
are calculated with respect to ^ . 



negative 5* results in the stau mass that increases as we run down from the GUT 
scale, allowing an acceptable stau mass even with mo = 0. 

The points B6 and B7 are representative of the f coannihilation region, and the 
breakdowns of their tunings are also shown in Table [5l The dependences on mo 
and mi/2 are similar to what was observed in the CMSSM. However, the dominant 
sensitivities are now to niH^^. For both these points the soft Higgs mass-squared 
terms are large and negative at the GUT scale. As we have seen, these soft param- 
eters have a significant effect on the stau RGE. Therefore the coannihilation strip 
exhibits tuning with respect to these parameters. The total sensitivity remains low 
suggesting that, even though the soft Higgs masses play a role in the running, the 
dominant contribution to the stau mass is still from Mi. 

More distinctive deviations from the familiar CMSSM phenomenology arise in the 



19 



forms of the strong vertical dark matter regions at particular values of mi/2- In 
panel (a) three vertical strips are present. To understand these lines we need to 
consider the mass and composition of the lightest neutralino. The bino component 
of the lightest neutralino is determined by Mi{EW) ~ 0.4mi/2(Gt/T), whereas the 
wino component is determined by M2{EW) ~ 0.8mi/2{GUT). Hence, M2{EW) > 
Mi{EW) throughout the NUHM parameter space, and we never have a large wino 
component in the LSP. Of more importance is the higgsino component, determined 
by fi{EW). When fi{EW) ~ Mi{EW), there will be a sizeable higgsino component 
in the LSP. In panel (a) we have set fi = 400 GeV and rriA = 400 GeV. Therefore, 
when mi/2 ~ 1000 GeV, Mi{EW) ~ fi and the lightest neutralino will be a 
bino/higgsino mixture. However, for mi/2 ^ 1000 GeV the lightest neutralino is 
mainly a higgsino, with a mass m^o ^ 400 GeV, whereas for mi/2 ^ 1000 GeV 
the Xi is predominantly a bino and has a mass determined by Mi{EW). 

With this in mind, we can understand the vertical lines in panel (a) at particular 
values of mi/2. At mi/2 = 500 GeV, the lightest neutralino is a bino with a mass 
m^o ^ 200 GeV. As the pseudoscalar Higgs mass is = 400 GeV throughout, 
this results in resonant neutralino annihilation through the pseudoscalar Higgs. 
As a result, the relic density is below the WMAP value across the region 450 < 
mi/2 < 530. On the edges of the resonance the relic density may fall within the 
narrow range favoured by astrophysics. The lines are mostly plotted in purple and 
blue, corresponding to a large degree of dark matter fine-tuning A*^ = 30 — 300. 
However, it is interesting that the edge of the resonance at larger mi/2 is plotted in 
red. This is the first instance of an acceptable pseudoscalar Higgs resonance region 
with relatively low dark-matter fine-tuning, thanks to the larger higgsino fraction 
in the LSP at larger mi/2. Both the LSP mass and the mass of rriA are sensitive 
to /i, so the mass of the LSP and rriA are coupled. This mitigates the dark matter 
tuning of the pseudoscalar Higgs funnel to an extent. 

Points Bl and B2 lie on either side of the pseudoscalar Higgs funnel, where it is 
interesting to break the dark matter fine-tuning measure down into its component 
parts. Unsurprisingly, the dark matter fine-tuning is due to a balancing act between 
the mass of the pseudoscalar Higgs, as shown by the large sensitivity to , and 
the neutralino mass, as shown by the sensitivity to mi/2- 

The other vertical strip is a wide region around mi/2 ~ 750 GeV, where the pre- 
dominantly bino LSP acquires a sufficient higgsino admixture to suppress the relic 
density to the observed value. This wide band is plotted in red and purple, repre- 
senting a tuning of 10 — 30, similar to that of the low-mo ^nd of the bino/higgsino 
region within the CMSSM. There is also a thin line of green at the base of this 
band where it meets the coannihilation strip. This shows that the interaction of a 
stau coannihilation region with a bino/higgsino LSP reduces the overall tuning of 
either region alone. 
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Point B3 is representative of this region. The primary dark matter fine-tunings are 
clearly with respect to m//^ and mi/2- In this case, the dark matter fine-tuning is 
related to the balance between the roles of these terms in determining fi and Mi at 
the electroweak scale. This is the first instance of a mixed bino-higgsino region at 
low mo that we find in the NUHM. It is therefore interesting to see that, even with 
a TeV-scale value of m^^j, the tuning remains relatively small. This is in contrast 
to the point A3, at which a TeV-scale soft Higgs mass gave large fine-tuning. 
To understand the origin of both electroweak and dark matter fine-tunings, it is 
useful to consider the analytic form of the dependence of the low-energy parameter 
fi{EW) on the GUT-scale soft inputs. For tan/5 = 10, we have: 

2 

^ = -0.9V + O.OlOm^^ - 0.19M2 - O.OOITM^^ - O.GSm^^ + 0.38m^3 

+0.38m^3 + 0.093^2 - O.OllAMi - O.O7OAM2 - O.3OAM3 
+2.51M| + O.OO59M1M2 + O.O28M1M3 + O.I95M2M3, (17) 

from which we can see that, when ttlq^ = rn'j^_^^, the terms from the scalars ap- 
proximately cancel. This explains the jump in sensitivity when this universality is 
broken. On the other hand, in order to obtain a small value of /i, these soft scalar 
terms need to provide a large contribution to balance out the contribution from 
M|. This explains why we find a higgsino/bino only at large mo in the CMSSM. 
In the NUHM, it is unnecessary to go to large mo, just large frijj^. By keeping 
mo, and thus ttlq,^, small one keeps the dark matter and electroweak fine-tunings 
associated with these parameters under control. The remaining large electroweak 
fine-tuning associated with point B3 is due to M1/2 being quite sizeable. 

In panel (b) we have /i = 400 GeV and m^ = 700 GeV. Here we see only the lower 
edge of the pseudoscalar resonance. This band lies at 750 GeV and is plotted in 
purple, once again showing the large dark matter fine-tuning we expect of such 
resonances. The upper edge of the resonance never appears, because at larger mi/2 
the lightest neutralino is dominantly higgsino. Therefore, at values of mi/2 above 
the resonance, the higgsino nature of the LSP suppresses the relic density so that 
the relic density never rises enough to fit the WMAP measurement. 

Point B4 illustrates the component dark matter fine-tunings at the resonance. As 
with point Bl, the dark matter fine-tuning is large, and primarily due to sensitiv- 
ities to m^^ and mi/2. These are due to their effects on the pseudoscalar Higgs 
mass and the mass of the LSP, respectively. 

In panel (c) we set /i = 700 GeV and mA = 400 GeV. As the pseudoscalar mass is 
the same, the pseudoscalar Higgs funnel is centred in the same place as in panel (a). 
However, the higgsino fraction of the LSP has dropped, increasing the overall relic 
density. Therefore one must go closer to the resonance before the enhancement to 
the annihilation cross section is sufficient to fit the observed relic density. Therefore 
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the WMAP lines are closer to the peak of the resonance, and require more precise 
dark matter fine-tuning than for the lower value of /i. This is shown by examining 
the component dark matter fine-tunings for point B5. The dark matter fine-tuning 
is large, and primarily due to m/^-j through its infiuence on the pseudoscalar Higgs 
mass. In panel (c) there is no region in which the LSP is higgsino, due to the larger 
value of fi. One would have to go to mi/2 ~ 1700 GeV before the LSP acquires a 
significant higgsino fraction as /i = 700 GeV. 

Finally, in panel (d) we take fi = 700 GeV and = 700 GeV, and we see the 
pseudoscalar Higgs resonance at mi/2 = 800 GeV as expected, and it remains finely 
tuned, as before. As in panel (c), the value of fi is too large to find a region of the 
parameter space in which the neutralino is higgsino. 

Overall, Fig. [3] shows some similarities and some marked deviations from the phe- 
nomenology of the CMSSM. A f coannihilation band appears in roughly the same 
region of the parameter space and exhibits slightly larger tuning. This increase is 
due to the effect of the soft Higgs masses squared on the f running. A further de- 
viation from the CMSSM comes in the form of a pseudoscalar Higgs funnel at low 
tan (3. In most cases this exhibits dark matter fine-tuning similar to the CMSSM, 
supporting the observation that resonances require significant dark matter fine- 
tuning wherever they appear. The exception is when the LSP has a significant 
higgsino fraction. The sensitivity of both and m^o to fi lowers the required 
dark matter fine-tuning significantly. There also appears a band of higgsino/bino 
dark matter and we find that the tuning is of the same order as in the CMSSM, 
which is surprising. In the CMSSM the focus-point region has relatively low dark 
matter fine-tuning because of the cancellation among the scalar masses. By break- 
ing this universality, this cancellation is broken. However, it also allows us to 
access mixed bino/higgsino regions at low rriQ. This enhances other annihilation 
channels, such as that via t-channel slepton exchange, and we require a smaller 
higgsino admixture to obtain a suitable dark matter relic density. These regions 
recover the relatively low dark matter fine-tuning of the focus-point region of the 
CMSSM. 

Finally, we note that we have also calculated the electroweak fine-tuning across 
these planes in the NUHM parameter space, and find it to be very similar to that 
of the CMSSM planes studied previously. 

We now consider the behaviours of these regions as tan /3 increases. In Fig. H] we 
take fi = 400 GeV, = 700 GeV, Aq = and sign(/i) positive, and in panels 
(a)-(d) we take successively larger values of tan/3[§. 

As the value of tan f] increases, the mass and composition of the neutralino across 
the (mi/2,?Tio) plane remains essentially unaltered. Therefore we find the lower 

6 Note the CMSSM point in panel (a) of Fig. H 
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Figure 4: Sample {mi/2,mQ) planes of the NUHM parameter space with Aq = 0, fi = 
400 GeV, mA = 700 GeV and sign(;u) and the following values of tan/3: (a) tan/3 = 10, 
(b) tan/3 = 20, (c) tan/3 = 35, (d) tan/3 = 50. The Roman cross in panel (a) shows the 
single point where the NUHM makes contact with the MSSM. 



edge of the pseudoscalar funnel in the same place in all panels. Studying the 
breakdown of the dark matter fine-tunings of the pseudoscalar Higgs funnel for 
large tan j3 at point B8, we find that in this case all parameters show a significant 
tuning. The primary tunings come from m^/j and tan/3. From (fT0l) .(l8l) we can see 
that this comes from the dominant sensitivity of m\, through fi, on both tan/5 and 

The pseudoscalar Higgs funnel remains in the same place as tan j3 increases, whereas 
the f-coannihilation strip moves considerably. For larger tan /5 the Yukawa contri- 
bution to the stau running is enhanced. This suppresses the stau mass, requiring 
larger mo to avoid a stau LSP. For very low mo there is a brick red region. This is 
not due to a failure of REWSB, but rather it is due to a tachyonic stau mass. 

By increasing tan (3 we reduce the dominance of Mi in the running of the stau mass- 
squared. This breaks the link between the masses of the f and the Xi^ resulting 
in an increase of the tuning required for the coannihilation region at larger values 
of tan/3. Point B9 is representative of the f coannihilation region at large tan/3. 
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Table 6: Points B8 and B9, shown in Fig. H] exemplify the behaviour of the pseudoscalar 
Higgs funnel (B8) and the stau-coannihilation region (B9) at large tan/3 within the 
NUHM. We present a breakdown of the dark matter fine-tuning with respect to each 
parameter of the NUHM. We give the value of m'jj^ ^ , although the tunings are calculated 
with respect to niHi 2- 

Though there is large sensitivity to the soft gaugino mass mi/2, this is eclipsed by 
the sensitivities to tan/3 and rriQ. 

4.4 Sample (/i, m^) planes 

We now discuss other planes so as to explore more features of the NUHM model. 
In Fig. Owe consider m^) planes with Aq = 0, tan/3 = 10 and sign(/x) positive, 
taking different values of mo and mi/2 in the various panels. As before, though we 
plot results in terms of fi and rriA, the primary variables are m|f^ ^. 

We first consider some overall features. Low fi is ruled out in all cases by the 
appearance of a light chargino (light blue). In panels (a) and (c), the low value 
of mi/2 = 300 GeV results in a Higgs boson with rrih < 111 GeV across much of 
the parameter space. This bound is very sensitive to variations in the top mass, 
so regions that are excluded here could be allowed with a higher top mass. On 
the other hand, in panels (b) and (d), only very low values of rriA result in a 
problematically light Higgs. In panels (a) and (b), the low value of niQ = 100 GeV 
results in a light f. This gives a region at low n and rriA in which the f is the 
LSP and as such is ruled out (light green). Finally, in panel (a) mi/2 and mo are 
light enough to give light sneutrinos. For large fi and m^, one of the sneutrinos 
becomes the LSP, ruling out this corner of the parameter space 0. 

^Sneutrinos are massive, neutral and weakly-interacting, and so could in principle account for 
the dark matter. However, they generally give VLcDuh^ ^ ^'cdm^ ■ This can be avoided in 
some models with right-handed neutrinos, but this possibility lies beyond the scope of the MSSM, 
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Figure 5: Sample NUHM {fi, hia) planes with Aq = 0, tan/3 = 10 and sign(^) positive, 
and different values of mo and 'mi/2- (a) rriQ = 100 GeV, mi/2 = 300 GeV, (b) niQ = 
100 GeV, mi/2 = 500 GeV, (c) mo = 300 GeV, = 300 GeV, (d) mo = 300 GeV, 

^1/2 = 500 GeV. The Roman crosses in each panel show where the NUHM meets the 
CMSSM. 

Before considering the dark matter regions, we first note the composition of the 
hghtest neutrahno in different regions of the parameter space. In panels (a) and 
(c) Mi{EW) ^ 120 GeV, whereas in panels (b) and (d) M^iEW) ^ 200 GeV. 
Therefore for fi ^ 120(200) GeV the LSP is predominantly a bino with a mass 
m^o ^ 120(200) GeV. Below these values of n, the lightest neutralino acquires a 
significant higgsino fraction. 

As before, the regions that fit the relic density favoured by WMAP are displayed 
with the corresponding dark matter fine-tuning colours. The vertical dark matter 
band at low /i in all panels features a mixed bino/higgsino dark matter particle, 
and is plotted in red and purple, corresponding to a dark matter fine-tuning ^ 
20 — 40. Points CI and C9 are representative of the higgsino/bino band, and their 
dark matter fine-tunings are broken down in Table [71 As one would expect, the 
fine-tuning is due to a balancing act between rriHi and mi/2, as these determine 

so we do not consider it further here. 
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Table 7: Points Cl-11, shown in Fig. El illustrate the behaviour of mixed bino-higgsino 
dark matter (Cl,5,9), sneutrino coannihilation (C3), the pseudoscalar Higgs funnel (C4), 
stau-coannihilation (C2,6,7,ll), and the pseudoscalar funnel region with a mixed bino- 
higgsino LSP (C8,10). We present breakdowns of the dark matter fine-tuning with 
respect to each parameter of the NUHM. We give the value of rriff^ ^, but the fine-tunings 
are calculated with respect to ijih^ ^ • 
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the composition of the hghtest neutrahno, and thus its annihilation rate. 

Point C5 also represents a bino/higgsino region. It is unusual in that, normally, the 
larger the higgsino component, the more efficient the annihilation. However, in this 
case, as /i is lowered from point CI, the higgsino fraction increases and yet the dark 
matter relic density rises again. Indeed, C5 is almost evenly split between higgsino 
and bino components and yet it fits the WMAP relic density measurement. This 
is due to the fact that, as /x drops, the mass of the neutralino is lowered. At point 
C5 the LSP becomes lighter than the Z, shutting off the annihilation channel 
XiXi ~^ This reduces the annihilation cross section enough that the relic 

density is acceptable once more. 

Another region that is easy to pick out is the pseudoscalar Higgs funnel along 
niA = 240(400) GeV for mi/2 = 300(500) GeV respectively. As before, these 
funnels require significant dark matter fine-tuning and as such are predominantly 
plotted in blue and black, showing a tuning A*^ > 100. Point C4 is a representative 
point whose dark matter fine-tuning breakdown we display in Table [71 In previous 
pseudoscalar Higgs regions the dark matter fine-tuning was due primarily to rriHi 2 
and mi/2. Here we find that the sensitivity to the Higgs masses has increased 
significantly. From ( ITOl) this is easy to understand. If we increase fi while keeping 
rriA the same we must carefully balance the large m|^^ ^ contributions to give the 
required m^. This careful balancing manifests as a steadily increasing sensitivity of 
niA to the Higgs soft masses as we increase /x. This translates to a large sensitivity 
of the pseudoscalar Higgs funnel. 

At the other end of the spectrum, there is a region of the pseudoscalar Higgs funnel 
at low /i with remarkably low tuning. This occurs when there is a significant 
higgsino fraction in the LSP, such as at points C8 and CIO. In this region, both 
rriA and the neutralino mass are sensitive to /i. This results in the mass of the 
neutralino and the pseudoscalar being coupled, and reduces the sensitivity of the 
mass difference = m^ — 2m^o. At points C8 and CIO the dominant annihilation 
channels are to heavy quarks via an s-channel pseudoscalar Higgs. Remarkably the 
total dark matter fine-tunings of the points are only 6.7 and 5.1 respectively. 

As the Ue^fj, become the LSPs in the large /i, large region of panel (a), there is a 
corresponding sneutrino coannihilation region lying parallel to its boundary, which 
is plotted in purple and blue indicating a dark matter fine-tuning A*^ > 80. Point 
C3 is a representative of this region, whose dark matter fine-tuning breakdown is 
also displayed in Table [71 The dark matter fine-tuning is large, and comes primarily 
from the Higgs sector. It is the existence of large negative that allows for light 
sneutrinos. Thus the sneutrino masses are very sensitive to the Higgs soft mass- 
squared parameters, and this is reflected in the dark matter fine-tuning. There 
is also some dark matter fine-tuning with respect to mi/2 that is typical of the 
need to balance the bino mass against that of a coannihilation partner with an 
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uncorrelated mass. 

Finally, the light f at low /i and mo has an effect on the dark matter relic density. 
As the mass of the f is reduced, the annihilation cross section is increased via 
t-channel slepton exchange. Also, as one approaches the region in which the stau 
is the LSP, there are additional contributions from f — x? coannihilation processes. 
These two effects combine to give dark matter bands along the edges of the stau 
LSP region in panels (a) and (b). Points C2 in panel (a) and C6,7 in panel (b) 
are representative points. At point C2 the annihilation proceeds through equal 
parts of t-channel e^, fiR, f annihilation (15-20% each), annihilation to b,b via 
off-shell pseudoscalar Higgs bosons (18%) and f coannihilation (15%). Only the 
coannihilation processes would be expected to exhibit a high sensitivity to the soft 
parameters, as t-channel processes are fairly insensitive and the point is far from 
the pseudoscalar resonance, reducing significantly the sensitivity of the s-channel 
pseudoscalar process. As a result, we have a region that arises from a mixture of 
channels and exhibits low tuning. The subdominant role of coannihilation explains 
why there is so little dark matter fine-tuning with respect to mi/2. The role of the 
stau in both the coannihilation and t-channel processes explains the dominant dark 
matter fine-tuning with respect to mo, and the dependence on mj^j appears from 
running effects. 

Unfortunately, point C2 results in a light Higgs with m/j = 110 GeV, which is 
probably unacceptably low, even allowing for the theoretical uncertainty in the 
calculation of its mass. On the other hand, panel (b) has a larger value of mi/2 
and hence Higgs mass. However, the masses of the LSP and the sleptons are 
also increased. This decreases the slepton t-channel annihilation cross sections, 
requiring larger contributions from processes that are finely tuned in order to fit 
the WMAP relic density, which is apparent at points C6 and C7|^. At point C6, t- 
channel slepton annihilation only accounts for 3% of the annihilation rate via each 
channel (9% overall). The remaining 91% is made up entirely of coannihilation 
processes, dominantly with f, but also cr, jxR. As this plane has low mo, the slepton 
masses are predominantly determined by mi/2- Once again, there is the familiar 
pattern of dark matter fine-tunings for a low-tan /3, low-mo slepton coannihilation 
region. The overall dark matter fine-tuning is low, and what fine-tuning does exists 
is due to mo and mi/2. Point C7 tells a slightly different story. The pattern of 
annihilation channels is almost identical to C6, and we see the typical dark matter 
fine-tunings of a coannihilation region in the sensitivity to mo and mi/2- However, 
the dark matter fine-tuning with respect to ttih-^^ has increased dramatically, due 
to the massive increase in rn\j^ ^ between points C6 and C7. Now the stau running 
is dominated by the Higgs mass-squared terms rather than the gaugino mass, and 
the coannihilation region becomes fine-tuned once again. 

^ We note that there is a CMSSM point very close to C6 
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There is one further interesting region. In panel (d) at low there is a kink 
in the higgsino/bino region. The band moves to larger and the dark matter 
fine-tuning drops dramatically. The band is plotted in green rather than purple, 
indicating a dark matter fine-tuning of less than 10. The kink in the band appears 
at niA = 280 GeV. Around this region the LSP is predominantly a bino with a small 
but significant higgsino component, and the LSP has a mass of around 200 GeV. As 
the pseudoscalar mass drops, the masses of the heavy Higgs, H, and the charged 
Higgses, H"^, also decrease. Around ijia = 280 GeV, the annihilation channels 
XiXi ~^ hAjW^H"^ , ZH open up, which are kinematically forbidden at larger 
rriA- These can progress through either t-channel neutralino (chargino) exchange 
or s-channel Higgs and Z processes. They require a small higgsino component, 
but significantly less than the higgsino/bino region represented by point C9. This 
balance of the higgsino and bino components of the LSP appears in the sensitivity 
of point Cll on mi/2 and m/^j. Thus Cll represents a higgsino/bino region with 
low dark matter fine-tuning - something that does not exist in the CMSSM. This 
is because a large negative mj^^ is needed to achieve low mA,H,H±- 

We now consider in Fig. [6] the behaviours of these regions as tan f] increases. We 
have set mo = 100 GeV, mi/2 = 300 GeV, Aq = and increase tan/3 in steps in 
each panel. 

We note first the bulk features of the plane. As noted previously, increasing tan P 
decreases the mass of the lightest stau. Thus plots at larger tan/3 have larger 
regions ruled out because the LSP is a f, and we do not show very large tan/3 
because at tan /3 = 50 the stau mass becomes tachyonic across the entire plane. 
By tan /? = 35 the light stau rules out all the parameter space below /i = 200 GeV. 
The mass of the light Higgs is also sensitive to tan (3, and is in all cases very close 
to rrih = 111 GeV, so it only takes a small shift to cause a significant change in the 
area plotted in light grey. 

The most significant change in the dark matter phenomenology is due to the varying 
f mass. Between panels (a) and (b) the stau bulk/coannihilation region moves to 
larger /i and rriA- We also find a significant stau region at large m^. These features 
are represented by points C15 and C14 respectively B Comparing C15 directly to 
C2, we see from Table [H] that the dark matter fine-tuning is due primarily to mo 
and tan/3. This is because these parameters determine the mass of the lighter 
stau and this is the primary source of sensitivity for bulk regions. There is also a 
degree of sensitivity to mi/2, as there is a significant coannihilation contribution 
that requires the LSP and stau mass to be balanced. At point C14 one has similar 
degrees of dark matter fine-tuning with respect to tan/3, mi/2 and mo. However, 
there is now also large dark matter fine-tuning with respect to rriHi 2 ? due to the 
larger magnitude of the soft higgsino mass-squared terms. The stau mass in this 

^ We note that there is a CMSSM point very close to C15. 
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Figure 6: Sample NUHM (/i, itia) planes with Aq = 0, tuq = 100 GeV, mi/2 = 300 GeV 
and sign(^) positive, and the following values of tan/3: (a) tan/? = 10, (b) tan/? = 20, 
(c) tan P = 35. We do not show a plot for tan /? = 50 as the parameter space is entirely 
excluded. The Roman crosses in each panel show where the NUHM meets the CMSSM. 

region becomes highly sensitive to m'jj^. 

The other regions are little changed from before. Point C12 exemplifies the mixed 
bino/higgsino region at increasing tan/3. It can be compared directly to point CI, 
and we see that the component dark matter fine-tunings are virtually identical. 
Point C13 is representative of the sneutrino coannihilation region and can be com- 
pared to point C3. Once again the dark matter fine-tuning is due primarily to the 
soft Higgs masses through their impacts on the running of the sneutrino masses. 

Much of the low-mo parameter space is forbidden by a light Higgs and/or a light 
stau. We now consider the effect of increasing tan/3 in a more open part of the 
parameter space. We take Fig. M^d) with mo = 300 GeV, mi/2 = 500 GeV as a 
starting point and increase tan j3 steadily, as seen in Fig. [71 In contrast to Fig. [6], 
the bulk features remain fairly stable for moderate values of tan/3. The first hint 
of a change appears in panel (c) at tan P = 35, where we see a small region at large 
rriA in which the stau is the LSP. This expands to cut off low fi for tan f3 = 50. 
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Table 8: Points C12-15, shown in Fig. [U are representative of the higgsino/bino region 
(C12), the sneutrino coannihilation region (C13) and the stau-coannihilation/bulk region 
(C15,14) with increasing tan/3 within the NUHM. We present a breakdown of the dark 
matter fine-tuning with respect to each parameter of the NUHM. We give the value of 
m^j^ 2 but the tunings are calculated with respect to rriHi 2 ■ 
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Table 9: Points C16-19, shown in Fig. [TJ illustrate the behaviours of the mixed bino- 
higgsino, the pseudoscalar Higgs funnel (C16) and the stau-coannihilation/bulk regions 
(C17,18,19) at increasing values of tan /3 within the NUHM. We present a breakdown of 
the dark matter fine-tuning with respect to each parameter of the NUHM. We give the 
value of m?ij^ ^ , but the tunings are calculated with respect to 2 • 



There are few dark matter surprises at larger tan (3. The pseudoscalar Higgs fun- 
nel and mixed higgsino/bino regions remain relatively unaltered throughout. The 
interaction of the pseudoscalar Higgs funnel with the higgsino/bino LSP continues 
to provide a favourable degree of tuning in panels (b) and (c). We take point C16 
as a representative point, and break the tuning down in Table As for point CIO, 
the tuning is small and the annihilation is primarily due to annihilation to heavy 
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Figure 7: Sample NUHM (/^, tua) planes with Aq = 0, ?no = 300 GeV, 7ni/2 = 500 GeV, 
sign(/i) positive and different values of tan/?: (a) tan/3 = 10, (b) tan/? = 20, (c) tan/? = 
35, (d) tan/3 = 50. The Roman crosses in each panel show where the NUHM meets the 
CMSSM. 

quarks via an s-channel pseudoscalar Higgs. 

Point C17 exemplifies the behaviour of a predominantly bino LSP with a small 
higgsino admixture that can annihilate to hA, ZH and W^H^. As with point 
Cll, the dark matter fine-tuning is small and mostly due to the composition of 
the LSP, through mi/2 and rriH^. 

Point C18 is in the new stau coannihilation region that appears at large tan (3. For 
mo = 300 GeV, mi/2 = 500 GeV the staus are too heavy to contribute significantly 
to t-channel slepton exchange, so this region is pure coannihilation. The stau 
mass is mainly determined by mo and tan/3, and must be balanced against a 
predominantly bino LSP. Therefore, the tuning is dominated by tan/3 and mo 
with a secondary dependence on mi/2- The coannihilation grows significantly by 
tan/3 = 50 and point C19 represents this trend. As with point C18, we find the 
tuning to be due to mo and tan/?, with a secondary dependence on mi/2- 

Throughout all of these parameter scans we have also calculated the electroweak 
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fine-tuning and found it to be of tlie same order as tliat found in tlie CMSSM for 
typical scales of soft masses considered. 



4.5 Sample (/i,mi/2) planes 

Finally,we consider sample (/i, mi/2) planes in the NUHM. These are interest- 
ing, e.g., because /x and mi/2 are the parameters that determine the mass and 
composition of the lightest neutralino E^. 




Figure 8: Sample NUHM (/x, mi/2) planes with Aq = 0, tan/3 = 10, sign(;u) positive 
and varying tuq and niA- (a) uiq = 100 GeV, = 500 GeV, (b) mo = 100 GeV, 
niA = 700 GeV, (c) mo = 300 GeV, = 500 GeV, (d) mo = 300 GeV, = 700 GeV. 
The Roman crosses in each panel show where the NUHM meets the CMSSM. 

In Fig. [H] we set Ao = 0, tan/3 = 10 and take discrete values of mo and m^. We 
see that either low n or low mi/2 results in a light chargino that violates particle 
searches (light blue). Low mi/2 also results in problems with a light Higgs (light 

^°Note that in the following plots mi/2 is the GUT-scale soft mass, whereas /x is the electroweak- 
scale Higgs term. This is in contrast to the plots of [T3] where the plots were in terms of M2 (EW) 
and n{EW). 
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grey with a black boundary). On the other hand, large mi/2 results in a neutralino 
with a mass above that of the stau (light green). The exception is low fi where the 
neutralino is a higgsino and m^o is insensitive to mi/2- In panels (a) and (b) we 
have mo = 100 GeV. This, combined with low mi/2 and large fi results in a region 
in which the LSP is a sneutrino (light green). 

We see once again the familiar dark matter regions of the previous plots. The 
pseudoscalar Higgs funnel appears as a pair of horizontal lines and exhibits large 
dark matter fine-tuning, and is characterized by the point D4 in Table [TUl Here we 
see that the large dark matter fine-tuning is due to the soft Higgs masses through 
their influence on rriA, and to mi/2 through its influence on the neutralino mass. 

The exception to this large dark matter fine-tuning is where the pseudoscalar funnel 
interacts with a higgsino/bino LSP and there is a small corner with low fine-tuning, 
as characterized by point D7. The annihilation here is mainly to heavy quarks via 
an s-channel pseudoscalar Higgs, and yet the total tuning is only 5.7. As noted 
previously, this relatively small dark matter fine-tuning comes from the common 
sensitivity of and m^o on /i. 

There is also a f coannihilation region in all four plots, which lies alongside the 
region ruled out due to a stau LSP. It exhibits similar tuning to the CMSSM. We 
break down the dark matter fine-tunings of this region at points D2 and D3, finding 
that at both points the tuning with respect to mo and mi/2 is standard for a stau 
coannihilation strip at low mo 0. Point D3 has larger tuning because this region 
of parameter space requires large negative soft Higgs masses, which now dominate 
the determination of the mass of the light stau. 

The sneutrino coannihilation region shows up alongside the sneutrino LSP region. 
Once again we find it to require significant dark matter fine-tuning, although this 
decreases steadily as one moves to lower /i. Point D5 is a representative point with, 
as before, large dark matter fine-tuning that depends on the soft Higgs masses. 

Each plot also has a dark matter region at low /i that lies along a diagonal 
in the (/i,mi/2) plane, incorporating points Dl,6,9. These regions are mixed 
bino/higgsino regions. In all cases the pseudoscalar Higgs and heavy Higgs bosons 
are sufficiently massive that annihilation of the mixed LSP proceeds mainly through 
the channels XiXi ~^ W~^W~{ZZ), via t-channel chargino (neutralino) exchange. 
This process is very sensitive to the composition of the LSP and the masses of the 
exchanged particles. Therefore there is significant dark matter fine-tuning with 
respect to mj^j and mi/2 at all these points. 

Finally, we consider the point D8 where the coannihilation strip and the mixed 
bino/higgsino strips meet. The combination of annihilation channels has a benefi- 
cial effect, with the overall dark matter fine-tuning dropping to 9. 

^^This is also true for the CMSSM point seen in panel (b). 
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Tablp 10: Points Dl-9, shown in Fig. [HI are representative of bino-higgsino dark matter 
(Dl,6,9), stau-coannihilation (D2,3,8), the pseudoscalar Higgs funnel (D4) and its inter- 
action with mixed bino-higgsino dark matter (D7), and sneutrino coannihilation (D5). 
We present breakdowns of the dark matter fine-tuning with respect to each parameter 
of the NUHM. We give the value of rrijj^ ^ but the tunings are calculated with respect to 
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Figure 9: Sample NUHM (/i, mi/2) planes with Aq = 0, mo = 100 GeV, niA = 500 GeV, 
sign(/i) positive and tan fj varying: (a) tan [i = 10, (b) tan (3 = 20, (c) tan (3 = 35. We do 
not show a plane for tan [3 = 50, as this part of the parameter space is entirely excluded. 
The Roman crosses in each panel show where the NUHM meets the CMSSM. 

Once again it is interesting to go beyond tan/3 = 10, to understand how the 
phenomenology changes with tan/5. In Fig. [9] we consider (/i, mi/2) planes with 
mo = 100 GeV, = 500 GeV and steadily increasing values of tan/5. As we 
saw before, increasing tan (3 decreases the f mass, causing the stau LSP regions to 
encroach on the parameter space. By tan /5 = 35 the light stau rules out all values 
of low /i. As noted earlier, at such a low value of mo, tan/5 = 50 has a tachyonic 
stau and so is not shown here. 

The change in the stau mass is the dominant factor that changes the dark matter 
phenomenology. With the lighter stau, the contribution to neutralino annihilation 
from t-channel stau exchange increases. We consider two points DIO and Dll in 
panel (b). At point DIO the annihilation is still dominated by coannihilation effects, 
but the growing contribution from t-channel stau exchange helps to lower the dark 
matter tuning. The dark matter fine-tuning is predominantly due to mo and tan (3 
through their influence on the mass of the lighter stau, with a subsidiary fine-tuning 
with respect to mi/2. In contrast, point Dll lies in a dark matter band where the 
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Table 11: Properties of points DlO-13, shown in Fig. [9] which are representative of the 
pseudoscalar Higgs funnel (D13) and the stau-coannihilation/bulk region (DIO, 11, 12) at 
increasing tan/3 within the NUHM. We present a breakdown of the dark matter fine- 
tuning with respect to each parameter of the NUHM. We give the value of m'j^_^ ^ but the 
tunings are calculated with respect to tuhi 2 • 



annihilation of neutralinos is dominantly through t-channel slepton exchange. As 
a result the dark matter fine-tuning is small, and due primarily to niQ and tan/3 
through their infiuence on the slepton masses. 

As we move to larger tan/3, the coannihilation and bulk regions meet. In panel 
(c) we take point D12 as a representative of the meeting of these two regions. 
However, by this stage one needs large soft Higgs mass-squared parameters and 
the stau mass is sensitive to these, rather than to mo and mi/2. Therefore there is 
large fine-tuning with respect to niHi- Finally, point D13 is representative point of 
the pseudoscalar Higgs funnel for large tan jS. As before, we find the dark matter 
fine-tuning to be large and predominantly due to the soft Higgs masses. This 
chimes with the general behaviour of the pseudoscalar Higgs funnel throughout 
our study. 

5 Conclusions 

We summarize in Table [12] the ranges of dark matter fine-tunings found in the 
different dark matter regions appearing in the NUHM, which may be compared to 
those found previously in the CMSSM, as shown in Tables UfH and summarised in 
Table El Comparing first the bulk regions, which require the smallest amounts of 
dark matter fine-tuning in both the CMSSM and the NUHM, we see that CMSSM 
point A2 has a low amount of dark matter fine-tuning that is at the end of the range 
found in the NUHM. However, point A2 has too small a value of m/j, and hence 
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Bino/higgsino region, m^o < mH,A 


2-10 


V — x\ coannihilation 


15-200 



Table 12: A summary of the different dark matter regions present within the NUHM 
and typical values of the corresponding dark matter fine-tunings. We also note that 
combinations from many channels decrease the overall tuning. 



is not acceptable, whereas the NUHM can circumvent this restriction. Thus, the 
NUHM provides access to dark matter which is less fine-tuned than in the CMSSM. 
Turning then to the f coannihilation regions, we see that the NUHM offers an 
option of lower fine-tuning than that found in the CMSSM at point A5. For 
reasons explained earlier in the text, there are very large variations in the amounts 
of dark matter fine-tuning required in the pseudoscalar Higgs funnel region of the 
NUHM, and the amount found at the CMSSM point A4 lies within this range. In 
addition, the NUHM contains a region in which a bino/higgsino LSP annihilates 
via a pseudoscalar Higgs boson that requires dramatically less fine-tuning than 
the pseudoscalar funnel in the CMSSM. Likewise, NUHM analogues of the focus 
points A3 and A6 shown in Table [2] require dark matter fine-tunings that are 
substantially less than in the CMSSM. One new region appears in the NUHM that 
has no CMSSM analogue, namely the v — x\ coannihilation region. At least in 
the examples studied, this requires rather more dark matter fine-tuning than the 
f — x\ coannihilation region. 

Although it was not the primary focus of this paper, we have also calculated the 
electroweak fine-tuning across the NUHM parameter space, and found it to be of 
the same order of magnitude as in other MSSM studies. 

Generally speaking, the fact that the NUHM has more parameters offers more pos- 
sibilities to find regions with particularly small (or large) dark matter fine-tunings. 
The smaller amounts of dark matter fine-tuning generally occur in regions where 
several different (co) annihilation processes contribute to the final dark matter den- 
sity, e.g., where a bino/higgsino band meets a stau coannihilation region. Con- 
versely, there are instances where a tight correlation is necessary between two a 
priori independent MSSM parameters, such as the stau and neutralino masses along 
a coannihilation strip, which is imposed automatically in the lower-dimensional 
parameter space of the CMSSM, resulting in smaller amount of dark matter fine- 
tuning than might otherwise have been expected. 



38 



There have been several studies of the imphcations of prospective LHC and/or ILC 
measurements for the accuracy with which the astrophysical dark matter density 
could be calculated on the basis of collider measurements pT| . These studies have 
emphasized relatively favourable points in the CMSSM coannihilation region at 
low mi/2, where the dark matter fine-tuning is relatively low and the prospective 
collider measurements relatively accurate. Our analysis offers some prognosis as to 
which NUHM regions might be favourable for extensions of these analyses. Clearly, 
the presence of NUHM points with low dark matter fine-tuning and relatively light 
sparticles is encouraging a priori. However, we note that in general the NUHM 
relic density is sensitive to the separate and independent values of the soft Higgs 
masses mHi2- Detailed consideration of their determinations using collider data 
has not yet been given, as far as we know, and we recall that the prospects for 
measuring directly the masses of the heavier Higgs bosons at the LHC and ILC 
are limited, though the prospects of the latter would be improved with the 7 — 7 
option, or by going to higher energies as at CLIC. 

The extension of these collider dark matter studies to the NUHM is a large task 
that lies beyond the scope of this paper. However, this exploratory study has 
revealed some of the prospective opportunities and pitfalls. 
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